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ABSTRACT 

The bimodal galaxy distribution in the optical colour-magnitude diagram (CMD) 
comprises a narrow "red sequence" populated mostly by early-type galaxies and a 
broad "blue cloud" dominated by star-forming systems. Although the optical CMD 
allows one to select red sequence objects, neither can it be used for galaxy classification 
without additional observational data such as spectra or high-resolution images, nor to 
identify blue galaxies at unknown redshifts. We show that adding the near ultraviolet 
colour (GALEX NUV A eff = 227 nm) to the optical (g-rvs M r ) CMD reveals a tight 
relation in the three-dimensional colour-colour-magnitudc space smoothly continuing 
from the "blue cloud" to the "red sequence". We found that 98 per cent of 225 000 
low-redshift (Z < 0.27) galaxies follow a smooth surface g — r = F(M r , NUV — r) 
with a standard deviation of 0.03-0.07 mag making it the tightest known galaxy 
photometric relation given the ~ 0.9 mag range of fc-corrcctcd g — r colours. Similar 
relations exist in other NUV-optical colours. There is a strong correlation between 
morphological types and integrated NUV — r colours of galaxies, while the connection 
with g — r is ambiguous. Rare galaxy classes such as E+A or tidally stripped systems 
become outliers that occupy distinct regions in the 3D parameter space. Using stellar 
population models for galaxies with different star formation histories, we show that 
(a) the (NUV — r, g — r) distribution at a given luminosity is formed by objects having 
constant and exponentially declining star formation rates with different characteristic 
timescalcs with the red sequence part consistent also with simple stellar population; (b) 
colour evolution for exponentially declining models goes along the relation suggesting 
a weak evolution of its shape up-to a redshift of 0.9; (c) galaxies with truncated 
star formation histories have very short transition phase offset from the relation thus 
explaining the rareness of E+A galaxies. This relation can be used as a powerful galaxy 
classification tool when morphology remains unresolved. Its mathematical consequence 
is the possibility of precise and simple redshift estimates from only three broad-band 
photometric points. We show that this simple approach being applied to SDSS and 
GALEX data works better than most existing photometric redshift techniques applied 
to multi-colour datasets. Therefore, the relation can be used as an efficient search 
technique for galaxies at intermediate redshifts (0.3 < Z < 0.8) using optical imaging 
surveys. 

Key words: galaxies: (classification, colours, luminosities, masses, radii, etc.) - galax- 
ies: photometry - galaxies: stellar content - galaxies: distances and redshifts 



1 INTRODUCTION 



*The best-fitting photometric relation coefficients and other 
supporting technical information arc available at the project web- 
site: |http:/ /specphot. sai.msu.ru/ galaxics / 

fE-mail: chilOsai.msu.ru (1(J); iz@sai.msu.ru (1Z) 



Understanding observational aspects of galaxy evolution re- 
quires to classify them in regard to various properties such as 
morphology, luminosity, stellar population characteristics, 
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internal dynamics. In the present-day era of deep wide-field 
imaging surveys, there is a need in efficient mechanisms of 
galaxy classification and selection using minimal available 
information. Colour-colour and colour-magnitude diagrams 

(CMD) have been traditionally used f or this purpose. 

In the optical CMD (g - r,MA dstrateva et al.l l200ll : 



opt 

iBaldrv et all |2004 iBlanton et at! l2003al) , the very narrow 
"red sequence" (jVisvanathan fc Sandageilf977h (a(g - r) « 
0.04 mag) formed mostly by elliptical and lenticular galaxies 
is used to identify early-type members of galaxy clusters 
because at low redshift it moves as a whole remaining tight 
in the colour space. However, optical CMDs cannot be used 
for detailed classification of galaxies, either for the selection 
of other than red galaxies because of several degeneracies: 
(I) there is no unambiguous connection between the galaxy 
morphological type and its position on the CMD; (2) the 
red part of the CMD is contaminated by ~ 25 per cent with 
late-type galaxies having weak ongoing star formation (SF) 
attenuated by the dust; (3) the blue clo ud overlaps with 
the lo ci of "E+A" poststarburst galaxies (|Dressler fc Gunnl 
If 9831 ) (PSG) having b lue colours, early-type m orphology, 
often disky kinematics (| Chilingarian et al]|2009br i but no or 
weak ongoing SF. 

In the (NUV — r,M r ) space, both the red sequence 
and the blue cloud become pro nounced but quite broad 
(a(NUV-r) ~ 2 mag) sequences ijWvder et al.ll2007f ). Such 
a width of the red sequence is due to the UV flux sensitiv- 
ity to even small fractions of young stars that was shown 
to be connected to the environment of early- type galaxies 
jKavirai et alj|2007ft . At the same time, (I) the sequences 
are too broad to use them for the efficient photometric selec- 
tion of galaxies; (2) there is still an ambiguity between the 
NUV — r colour and a galaxy morphological class as well 
as the presence of ongoing SF: PSGs still reside in the blue 
cloud. 



2 THE UV-OPTICAL GALAXY 
PHOTOMETRIC SAMPLE 

2.1 Catalogue construction 

Using Virtual Observatory data mining, we constructed a 
photometric sample of ~225 000 galaxies excluding quasars 
and bright active galactic nuclei (AGN) based on their spec- 
tral clas sification by the Sloan Digital Sky Survey Data Re- 
lease 7 l|Abazaiian et al1l2009l ) in the absolute magnitude 
range -25 < M z < -15 mag at low redshifts (0.007 < Z < 
0.27). We cross-identified the spectral sample of SDSS DR7 
galaxies with the U V Galaxy Evolution Explorer satellite 
(Martin et alll2005l) Release 5 ( G ALEX GR5) catalogue in 
the CAS Jobs jSzalav et ai]l2002f ) catalogue access systems 
of SDSS and GALEX and rejected the matches separated 
by more than 3 arcsec on the sky. 

First, we employed the SDSS DR7 CasJobs service to 
select galaxies in the redshift range from 0.007 to 0.27 from 
the SDSS DR7 spectroscopic sample in the stripes covered 
(already or in the survey plan) by the Unit ed Kingdom In- 
frare d Telescope Deep Imaging Sky Survey (|Lawrence et al.l 
120071 ) Large Area Survey Data Release 8 (UKIDSS LAS 
DR8). We selected only the objects classified as galaxies by 
the SDSS spectroscopic pipeline (SpecClass — 2), that al- 



lowed us to reject quasars and prominent broad-line active 
galactic nuclei. This list contains 377 923 sources. 

After that we appended both GALEX GR5 and 
UKIDSS DR8 data by joining our reference list of objects 
with these surveys using the spatial match criterion, namely 
the best match within the angular separation of 3 arcsec. For 
the GALEX data we made use of the pre-c alculated cross- 
matc h between SDSS and GALEX surveys l|Budavari et al.l 
120091 ) accessible through GALEX CasJobs as the xsdssdr7 
table restricting the angular separation to < 3 arcsec. For 
the UKIDSS LAS we queried multi-cone search program- 
matic access interface with effectively the same parameters 
for every object from our initial list. The SDSS-GALEX 
join returned 223 646 galaxies detected in NUV, 144 639 in 
FUV (A eff = 155 nm), and 136 781 in both filters. The 
SDSS-UKIDSS match contains 176 868 galaxies detected 
in the Y band, 178 806 in J, 187 789 in H, 188 221 in K, 
among them 158 578 in all four NIR bands. For 96 939 
of those galaxies we had photometric data from GALEX 
NUV, including 59 994 with GALEX FUV measurements. 
All the technical op erations on t ables were performed with 
the stilts software |Tavlor]|2006l ). 

We used SDSS Petrosian magnitudes (PetroMag_*), 
GALEX extended source calibrated magnitudes (nuvjmag 
and fuv_mag), and UKIDSS Petrosian magnitudes (Petro- 
Mag_*) to construct multi-wavelength spectral energy dis- 
tributions (SED). Here we notice, that even though SDSS 
model magnitudes (modelMag_*) generally have lower for- 
mally computed statistical uncertainties than Petrosian 
magnitudes, especially in blue photometric bands, in case 
of "blue cloud" galaxies they are often hampered by the dif- 
ferences between the observed light distribution and those 
assumed (axisymmetrical exponential or de Vaucouleurs) 
for the computation of model magnitudes. Petrosian mag- 
nitudes may underestimate the total galaxy flux by 15- 
20 per cent in ca se of face-on de Vaucouleurs profiles 
|Yasuda et al.ll200ll ). However, in our case this offset is simi- 
lar for SDSS and UKIDSS data while for the GALEX mea- 
surements it is not important because of high photomet- 
ric uncertainties significantly exceeding 15 per cent for red 
galaxies having the light profile shape affected by this effect. 
As long as GALEX and SDSS contain photometric mea- 
surements in the AB system, but UKIDSS magnitudes are 
in the Vega system, we ap plied zero-point tra nsformations 
available in the literature (|Hewett et al.l [20061 ) to the NIR 
magnitudes. We are using integrated photometry of galax- 
ies, therefore aperture effects have little importance in the 
present study and we do not need to apply aperture correc- 
tions. 

Then, all magnitudes were corrected for the effects of 
Galactic extinction. The UKIDSS and SDSS catalogues pro- 
vide selective extinction values in all photometric bands, 
while for GALEX w e used the provided E(B — V) value 
l|Schlegel et all 19981 ) and computed extinctions in UV bands 
assuming A NUV = 8.87- E(B - V), A FUV = 8.29-E(B-V). 

At this point we created a calibrated photometric sam- 
ple of low-redshift galaxies in 11 bands, from far-UV to NIR. 
It is easily reproducible at any workstation with the Inter- 
net access, however, due to the data access policy, the DR4 
latest public release of UKIDSS catalogues has to be used 
instead of DR8. 

Systematic uncertainties of SDSS point source photom- 
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etry do not exceed 1 per cent (|lvezic et alj 12004 ) in gr, 
whereas for extended sources they may be a few times larger. 
However, since the red sequence in the optical CMD of our 
sample constructed from galaxies populating a large area 
on the sky is as tight as 0.03 mag, we conclude that ei- 
ther the systematic errors on optical magnitudes are within 
this range, or they are strongly correlated between g and r 
bands so that they cannot hamper the results of our analy- 
sis. Statistical uncertainties of SDSS photometric measure- 
ments are generally better than 0.015 mag reflecting the 
spectroscopic target selection of SDSS: galaxies from the 
spectral sample are at least a few magnitudes brighter than 
the limiting magnitude of the photometric survey. At the 
same time, the median value of NUV magnitude uncertain- 
ties is as large as 0.15 mag across the whole sample. However, 
the range of NUV — r galaxy colours is ~ 7.5 mag compared 
to ~ 0.9 mag in g — r. Therefore, the relative "resolution" of 
our analysis per colour range is very similar in both colours. 



2.2 Computation of fc-corrections 

As we compare photometric measurements for galaxies at 
different redshifts, we have to correct them for the changes of 
effective rest-fram e wavelengths of filter bandpasses known 
as fc-corrections dOke fc Sandagel 1 19681 ; iHogg et al.l |2002| ; 
iBlanton fc RoweisboOTT T Their computation is an important 
step for obtaining the fully calibrated homogeneous dataset. 
Here we provide some details regarding the fc-correction 
computation in GALEX UV bands, while the procedure for 
optical and NIR fi l ters w as exhaustively described earlier 
(|Chilingarian et all 120101) . The importance of accurate fc- 
correction computation is illustrated by the fact that earlier 
studies of galaxies in t he (NUV — r,g — r) colour-colour di- 
agram (|Yi et alj|2005h did not report a sequence of galaxies 
which would be obvious if one took a galaxy sample in a 
narrow redshift range. 

Due to high sensitivity of UV fluxes to the recent SF 
and mass fractions of young stars as little as 1 per cent, the 
UV-to-NIR SED of a galaxy usually cannot be precisely rep- 
resented by a single simple stellar population, that is, a pop- 
ulation of stars of the same age and metallicity. Therefore, fc- 
corrections cannot be computed by the s ingle SSP fitting as 
it can be done in optical and NIR bands (|Chihngarian et aU 
2010). Other effects, such as a non-thermal emission from 
a moderate-luminosity active galaxy nucleus, or emission in 
certain spectral lines, create additional difficulties in the UV 
bandpasses. Most of these effects were tackled and success- 
fully taken into account using the nonnegative matrix fac- 
torization |BIanton fc Roweia 20071 ) . In the same work, the 



authors computed 5 synthetic template spectra representa- 
tive of different galaxies and galaxy components. 

Here we used the 2-step process to compute fc- 
corrections for our galaxies. First, we constructed a sub- 
sample including some 25 000 galaxies detected in all 11 
bands with high signal-to-noise ratios in the UV bands. We 
then fitted their SEDs using a non -negative linear combina - 
tion of 5 representative templates l|Blanton fc Roweis 2007) 
atten uated using the Cardelli extinction law 1 Cardelli et al.l 
1 19891 ) leaving the colour excess E(B — V) a free parameter. 



S econd, we use d the TOPCAlQ table manipulation soft- 
ware (|Tavloij [20051 ) to visualise obtained fc-corrections as 
functions of redshifts and variou s observed colours, as w e 
did for optical and NIR bands (|Chilingarian et all lioloh . 
Similarly, we found that UV fc-corrections could be pre- 
cisely approximated by low-order polynomial functions of 
redshifts and certain colours. The best filter combinations 
are NUV - g and FUV - u for the NUV and FUV bands 
respectively with standard deviations of the surface fitting 
residuals of about 0.08 mag and 0.15 mag. Then, we used 
these approximations to compute fc-corrections for all galax- 
ies in our sample. The newly obtained approximations of 
fc-corrections are available from the new version of the "fc- 
corrections calculator" servic^f]. 



3 THE COLOUR-COLOUR-MAGNITUDE 
RELATION AND ITS PROPERTIES 

We inspected the combined GALEX-SDSS dataset visually 
using TOPCAT in three dimensions (M z , NUV — r, g — r) 
and detected a very thin continuous distribution of both, 
blue and red galaxies along a smooth surface with very few 
outliers. Then we approximated it with a low order two- 
dimensional polynomial (see Appendix B for details). Given 
the range of observed g — r colour of ~0.9 mag, the dispersion 
of the g — r residuals that decreases from 0.07 to 0.03 mag 
going from blue to red NUV — r colours without significant 
dependence on the luminosity at M z < —17.5 makes it the 
tightest known photometric relation of galaxies. 

At lower luminosities, the dispersion of the residuals in- 
creases. In our case this can be explained by significantly 
lower number of objects due to the spectroscopic target se- 
lection algorithm used in SDSS and also by the poor quality 
of photometric measurements, because dwarf galaxies have 
lower mean surface brightness values than giants and, conse- 
quently, their magnitudes cannot be precisely measured in 
relatively shallow wide-field surveys that we used to con- 
struct our catalogue. An additional factor increasing the 
scatter is the peculiar motions of galaxies inside clusters and 
groups not taken into account which hamper the Hubble law 
distance estimates. 

The 3D distribution of galaxies and its projection onto 
the (NUV — r, g — r) plane are shown in Fig.[T][2] The choice 
of the z band for the M axis is not that important: the 
relation behaves similarly when using riz or NIR absolute 
magnitudes. We selected the SDSS z band for presentation 
purposes as the 2 luminosity range is slightly higher than in 
the r band. 

lYi et al.1 (|2005l ) presented a (NUV - r, g - r) colour- 
colour plot for galaxies and stars in their fig. 1. However, 
the galaxy photometric measurements were not properly fc- 
corrected as authors did not possess multi-band photometry, 
therefore no tight colour relation was revealed. 

Rather tight relation in the (NUV — r, g — r) colour- 
colour diagram fo r norm al gala xies was m e ntion ed in 
Schiminovich et al l (|2007l ). Then ISalim et all (|2005l ) and 
Haines et al.l ( 20081 ) attempted to use both, optical and near- 
UV information to study the star formation histories (SFHs) 



http : //www. star .br is . ac.uk/-mbt/topcat/ 
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Figure 1. The projection of the colour-colour-magnitude rela- 
tion on-to the colour— colour plane. The upper panel demonstrates 
the logarithm of the number density plot in grayscale with solid 
lines showing the relations for galaxies of constant luminosity de- 
rived from the best-fitting polynomial surface equation. Four bot- 
tom plots show fitting residuals in different magnitude ranges as 
density plots with dashed lines indicating their ±lcr levels, which 
are normalised to the maximum value in every NUV — r bin. 
Residuals for the two low-luminosity intervals are computed with 
coarser binning compared to the brighter galaxies in order to ac- 
count for lower object counts at those luminosities. Red sequence, 
blue cloud, and the loci of certain types of outliers are identified. 
The direction of internal extinction is shown with a vector. 



and dust effects and the effects of environment on the evo- 
lution of galaxies. However, the residual scatter of objects 
from this relation still remains high (an order of 0.15 mag 
in g — r) due to the dependence of both galaxy colours on 
luminosity. Adding the absolute magnitude as the third di- 
mension decreases the scatter by a factor of 3 in the absolute 
magnitude range (—25 < M r < —15 mag). It is remarkable, 
that the relation in the 3D space is followed by star-forming 
as well as passively evolving galaxies. 




Figure 2. The 3-dimensional distribution of galaxies in the 
(NUV — r, M z , g — r) space. The 3D plot presents the density 
distribution of 225 000 galaxies in the colour-colour-magnitudc 
space with the increasing density going from yellow to red, the 
best-fitting polynomial surface as a mesh grid immersed in it, and 
standard deviations of the fitting residuals shown as bars in the 
(NUV — r, M z ) plane with their colours cor responding to the 
magnitude ranges in Fig. [l] PSGs (Goto 2007), and compact el- 
liptical galaxies (cE) are shown as tetrahedra and cubes. The top 
face of the plot demonstrates their projected positions on-to the 
(NUV - r, M z ) plane. 



3.1 Connection to morphology 

To explore the connection of the colour-colour-magnitudc 
relation to galaxy morphology, we employed the morp holog- 
ical catalogue of SDSS galaxies (|Fukugita et alj|2007l ) avail- 
able through the Vizier servic^E which contains morphologi- 
cal types for 1 465 intermediate luminosity and giant galaxies 
at Z > 0.03 from our sample. We see a continuous change 
of morphological types along the surface in the NUV — r 
direction with a typical dispersion of 0.7. . .0.8 Hubble type. 
At the same time, the optical g — r colour turns to be a 
very bad morphological indicator: the red sequence region 
contains galaxies of all morphologies from ellipticals to Sc 
late-type spirals. 

The two-dimensional histogram of morphologies vs 
NUV — r colours is displayed in Fig. [3] One can see that the 
SOa and Sa galaxy types span a very broad range of colours 
demonstrating the difficulties of the visual classification of 
early-type disc galaxies. It corresponds to a simple linear 
correlation between the NUV — r and a Hubble type which 
can be expressed as: 

TYPE = 6.6 - 1.1 • {NUV - r), 

where the "TYPE" values will correspond to the Hubble 
types as < for E, 1 for SO, 2 for Sa, 3 for Sb, 4 for Sc, 
higher values for Irr. Red outliers (NUV — r > 4 mag) 
above the surface (yellow and red points in Fig. [3j mostly 
have later types, i.e. spiral galaxies, while outliers below it 
(blue and violet points in Fig. [3} have earlier types compared 
to galaxies on the sequence with similar NUV — r colours. 



http : //vizier . u- strasbg . f r / 
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Figure 3. Connection betwee n visually determined galaxy mor- 
phology jFukugita et al,|2007l ) and a NUV — r colour. Black con- 
tours correspond to galaxies following the relation (numbers are 
for the counts), while individual outliers beyond lcr are shown as 
crosses with the colours representative of the deviation from the 
relation in the g — r colour: yellow to red for galaxies above the 
relation, blue to violet for objects below it. 

The connection between the morphology and the lumi- 
nosity suggesting that more luminous galaxies have earlier 
morphological types is much looser and may be affected by 
the selection effects in our sample. 

3.2 Outliers from the relation 

We identify several classes of outliers from the relation com- 
prising ~2 per cent of the total sample (see Fig. [T][2]). 

(i) Early-type PSGs sele cted from the catalogue of H<5- 
strong galaxies (|Gotoll2007l ) (359 matches with our sample) 
populate a region ^0.15 mag below the surface in g — r 
spanning 3 < NUV — r < 5 mag colours explaining the 
nature of "blue early-type outliers" from the morphology- 
(NUV — r) relation described above. These are galaxies with 
truncated or multi-modal SFH where the last strong star for- 
mation episode has just been finished. The passively evolv- 
ing newly formed stellar population reddens much faster in 
the NUV — r colour than in the g — r one so that a PSG 
at first departs from the blue part of the relation (right in 
Fig.[TJ and notably later (after 2.5-3 Gyr) moves up increas- 
ing g — r towards the locus of red sequence galaxies. 

(ii) Compact e lliptical galaxies (|Chilingarian et al1l2009al ; 
iPrice et al.ll2009h are residing above the red sequence region 
of the colour-colour-magnitude relation at the low luminos- 
ity part. A few examples of new cE galaxies are shown in 
Fig El However, their colours are never redder than those of 
the most massive galaxies at the bright red sequence end. 
This fact is explained by their formation via severe tidal 
stripping of more massive progenitors, most likely early-type 
disc galaxies, by massive elliptical or cluster/group domi- 
nant galaxies. Progenitors of cEs are stripped in the inner- 
most regions of galaxy clusters, when the SF is ceased be- 
cause their interstellar medium is already removed by the 
ram pressure stripping created by the hot intergalactic gas 



(IGunn fc Gottll 19721 ). Depending on the previous SFH, these 
objects must reside either on the colour-colour-magnitude 
relation, or slightly below it, in the PSG locus. During the 
relatively fast tidal stripp ing process lasting about 1 Gyr 
(|Chilingarian et al]|2009al ). their stellar population proper- 
ties and colours change insignificantly while the mass and, 
consequently, the luminosity may decrease by a factor of 10 
or more, hence moving a galaxy off the relation if it was sit- 
ting on it. Thus, red cE colours are explained by high stellar 
metallicities inherited from their progenit ors which is con- 
firmed by detailed studies of nearby cEs l|Rose et al.l 120051 ; 
IChilingarian fc BergondteulCh . There may be some very rare 
intermediate-age cEs originating from PSGs whose colours 
will be bluer, however their passive evolution will quickly 
move them above the colour-colour-magnitude relation. 

(iii) Dusty starforming galaxies such as edge-on spirals 
are sometimes found above the flattened red part of the 
relation (NUV — r ;> 4 mag) also being consistent with the 
locus of late-type morphological outliers. Their positions are 
explained by the extinction vector direction shown in Fig. [T] 
If internal extinction is very strong then a galaxy is moved 
up-right in the diagram and may end up above the locus of 
red sequence galaxies. 

(iv) Galaxies with strong ongoing SF but yet small mass 
fractions of newly formed stars including ongoing and recent 
mergers may have very peculiar colours because of strong 
nebular emission lines and/or large quantities of dust. 

(v) Narrow-line AGNs having low contribution of their 
nuclei to the total light in the optical band and hence clas- 
sified as normal galaxies by the SDSS pipeline may have 
strong UV excess. We did not apply any particular filtering 
to our data to exclude these objects, therefore our sample 
may be slightly contaminated by them at a sub-per cent 
level. 

(vi) "Non-physical" objects: galaxies casually overlapping 
with either foreground stars or galaxies at different redshifts 
create outliers which may be located in almost any part of 
the parameter space except the region very red in (NUV — r) 
and very blue in (g — r). 



4 DISCUSSION 

4.1 Effects of stellar population evolution and 
internal extinction 

The distribution of galaxies in the colour-colour-magnitude 
space is governed by three factors: (1) stellar mass, (2) star 
formation and chemical enrichment histories including the 
ongoing SF, and (3) internal extinction. Therefore, there 
must be a connection between galaxy positions on the di- 
agram and their stellar population properties. We fitted a 
sub-sample of 133 000 SDSS DR7 spectra with simple stel- 
lar population (SSP) mod els using the NBURSTS technique 
l|Chilingarian et al. I l2007bl la1) and hence obtained their SSP- 
equivalent ages and metallicities. The NBURSTS technique 
includes the multiplicative polynomial continuum that ab- 
sorbs flux calibration errors and makes the fitting insensi- 
tive to the internal extinction in a galaxy. We also notice 
that SDSS spectra obtained in 3 arcsec wide circular aper- 
tures may not be representative of entire galaxies in case 
of strongly extended objects with notable gradients of the 
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Figure 4. Connection between galaxy positions on the colour-colour— magnitude relation and their SSP-cquivalent stellar populations. 
The slice of the relation in a narrow luminosity range (—23 < M z < —21.5 mag) is displayed in the colour-colour projection. Mean ages 
of stellar populations of 40,000 galaxies obtained from the fitting of their SDSS DR7 optical spectra are colour— coded. The evolutionary 
tracks of stellar population models without internal extinction for the solar metallicity and various SFH are ovcrplotted (see text). The 
colours of the tracks correspond to the ages of synthetic galaxies formed 12 Gyr ago, the ticks on the tracks are given every 1 Gyr. The 
NUV — r colours in the models are empirically corrected by —0.7 mag (see the text). 



stellar population properties. However, we can draw some 
qualitative conclusions. 

Even in the over-simplified case of SSP-equivalent pa- 
rameters, we observe a strong connection between aver- 
age stellar population properties of galaxies and their po- 
sition in the (M z , NUV — r,g — r) parameter space. Galax- 
ies sitting close to the best-fitting surface exhibit moder- 
ate metallicity gradient as a function of luminosity and al- 
most no variations in the colour-colour plane except the 
blue end of the sequence (NUV — r < 2.5; g — r < 0.5) 
where the metallicity quickly decreases. The luminosity- 
metallicity relation of early-type galaxies known to be re- 
spon sible for the tilt of their op tical colour-magnitude rela- 
tion jKodama fc Arimoto|[l997l ) similarly causes the tilt of 
the colour-colour-magnitude surface in its red part. 

The observed age effects in the colour-colour projec- 
tion are more important. At intermediate and high lumi- 
nosities, the age smoothly increases along the sequence from 
~500 Myr at NUV - r « 1.5, g - r w 0.3 to 13 Gyr 



at the red end. For low luminosity galaxies, the oldest 
SSP-equivalent ages of red galaxies decrease to 10 Gyr at 
M z = —18 mag being in accordance with the known anti- 
correlation of mean stellar population a ges with luminosities 



of dwarf elliptical galaxies in clusters (|van Zee et al J 12004 



Michielsen et al . 2008; Ch ilingarian et al.ll2008l ; IChilingarianl 
20091 ; ISmith et al.ll2009al lbl). Because of the same effect, the 



upper "edge" of the broad red sequence in the (NUV — 
r;M z ) CMD is strongly tilted at M z > -19 mag. 

At all luminosities, there is a notable age gradient across 
the sequence, that is, higher g—r colours correspond to older 
populations. Also, the dispersion of age estimates (logt) in- 
creases while moving towards blue colours with values to- 
tally uncorrelated with colours at the very blue end of the 
sequence corresponding to mean stellar ages t <C 500 Myr. 

In Fig. [4] we present a luminosity slice (—23.0 < M z < 
—21.5 mag) of the sample having the median SSP-equivalent 
metallicity [Fe/H] = -0.02 dex (([Fe/H]) = 0.13 dex) and 
the median redshift Z = 0.11 ((Z) — 0.03) corresponding to 
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the light travel time of ~ 1.4 Gyr. At this redshift range, the 
3 arcsec wide apertures enclose a large fraction of light from 
galaxies. The qualitative behaviour of mean stellar popula- 
tion ages at other luminosities is similar. 

We also present the evolutionary tracks for galaxies 
having various SFH types. The models were const ructed 
from the PEGASE.2 (|Fioc fc Rocca-Volmerangelll997l ') mod- 
els comput ed using the syntheti c low-resolution BaSeL stel- 
lar library l|Leieune et al.|[l997l ) for NUV — r colours, and 
g — r colours predicted by another family of stellar popula- 
tion models jVazdekis et aljfeoiol ) computed using a large 
Medium-resolution Isaa c Newton Telescope library o f em- 
pirical spectra (MILES. ISanchez-BIazquez et alj[2006h ■ The 
combination of t he two families o f stella r population models 
was essential, as iMaraston et alj (|2009l ) demonstrated that 
the offsets between predicted and observed colours of red 
galaxies in the SDSS photometric system was due to the 
nature of synthetic stellar spectra used to construct stellar 
population models. The proposed solution was to use models 
based on empirical stellar spectra. The MILES stellar library 
used to construct models presented here, has the best cov- 
erage of the stellar atmosphere parameters compared to all 
other existing published sources except the high-resolution 
ELODIE library which, however, has too narrow wavelength 
coverage making the computation of g and r colours impos- 
sible. The NUV — r colours for all models displayed in Fig. 0] 
were empirically corrected by —0.7 mag. Th is offset is prob- 
ably o f the same nature as that described bv lMaraston et all 
(III), however no models based on empirical stellar spectra 
are available in the NUV yet. 

The tracks shown in Fig. [4] were computed as follows. 
First, we computed colours and luminosities of simple stel- 
lar populations having [Fe/H] =0 dex using the PEGASE.2 
code for the ages of 30, 50, 100 Myr and further till 17 Gyr 
with a step of 50 Myr. Second, we computed g — r colours 
from the MILES-based models and interpolated them to the 
same age grid. The resulting SSP track is shown in Fig. Sl- 
its young part is strongly offset from the observed distribu- 
tion of galaxies towards red NUV — r then joining the main 
red sequence concentration at ages t > 5 Gyr. Third, we 
integrated the computed SSP luminosities in different pho- 
tometric bands up-to 12 Gyr using two families of an SFH: 

(a) exponentially declining star formation rate (SFR) with 
the three characteristic timescales texp = 1,2,4 Gyr; and 

(b) truncated SFHs: constant SFR until a given moment of 
time from the galaxy formation epoch (t^ T = 4, 6, 10 Gyr) 
followed by the immediate SF cessation and passive evolu- 
tion afterwards. 

Our galaxy evolutionary tracks are simpler than real 
galaxies because they do not include the intrinsic metal- 
licity evolution and other processes such as gas infall from 
filaments or satellites, mergers, etc., however, we can use 
them for some qualitative conclusions: 

(i) SSP models at low and intermediate ages (t < 5 Gyr) 
to not have any corresponding galaxies observed that sug- 
gests that (obviously) none of the massive galaxies in our 
sample was formed recently and quickly. However, it matches 
quite well the locus of the oldest red sequence galaxies. 

(ii) The main colour sequence can be explained by galax- 
ies formed immediately after the Big Bang (about 12 Gyr 
taking into account the median redshift of our sample) and 



having various types of a SFH. In its main part, the slope of 
the colour-colour relation and the direction of the internal 
extinction match each other very well. 

(iii) A galaxy having a constant SFR will end up near 
the low blue end of the sequence and can be moved up-right 
along it by the internal extinction. A family of exponentially 
declining SFH with different characteristic timescales form a 
curved sequence well corresponding to the observed relation. 
The NUV —r colour evolution pace at 3 < NUV — r < 5 mag 
anticorrelates with the SFR characteristic timescale (see also 
IWvder et all (|200?f ) for a similar plot in the (NUV — r, u — r) 
colour space without any selection on the luminosity). Mod- 
els for lower metallicities well reproduce the colour-colour 
relation at lower luminosities suggesting the universality of 
exponentially declining SFHs. We stress that this SFH type 
is not the only one that is able to explain the observed 
galaxy distribution in the colour-colour-magnitude space, 
however, it is the simplest model with the smallest num- 
ber of free parameters compared to other alternatives (e.g. 
multiple starbursts, an exponentially declined law with an 
additional burst). 

(iv) Galaxies having truncated SFHs have a very short 
transition phase on their way to the red sequence region last- 
ing about 1 Gyr after the SF cessation when their NUV — r 
colour reddens radically, by ~ 3 mag while the g — r change 
remains about 0.2 mag. The PSG locus below the main 
colour sequence is well matched by this transition phase and 
their rareness is consistent with a short duration of the tran- 
sition. 

(v) Dusty star-forming galaxies above the sequence at red 
NUV — r colours are also explained: they are moved up-right 
from the sequence following the direction of the extinction 
vector. 

(vi) The shapes of evolutionary tracks for galaxies with 
exponentially declining SFH clearly shows that the evolu- 
tion of a majority of galaxies goes along the relation during 
6-8 Gyr. Therefore, we would expect a weak evolution of the 
presented colour-colour-magnitude relation shape at least 
up-to a redshift Z ~ 0.9, although the distribution of galax- 
ies on it will evolve. This suggests it to be a unique search 
instrument for distant galaxy clusters using broad band gij 
images. 

4.2 Colour— colour— magnitude relations in other 
colour pairs 

We found similar photometric relations in other colour pairs, 
however they are more strongly affected by observational 
biases and galaxy evolutionary phenomena. Colour-colour- 
magnitude relations involving only optical colours are very 
tight because of strong degeneracies between the colours but 
for the same reason have very limited astrophysical applica- 
tions. For example, in the (u — r, g — r, M r ) colour space the 
u — r is close to ~ 2 X (g — r) at all luminosities, i.e. they 
are linearly dependent for most galaxies, hence virtually no 
information is added by the third dimension. 

4-2.1 Other near-UV-optical colour combinations 

The colour-colour-magnitude relation remains in place 
when other optical colour combinations together with the 
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NUV are used provided that there is enough wavelength 
lever in the optical colour to distinguish between red and 
blue galaxies. That is, colours like u — r,g — i, g — z,r — z 
but not r — i and i — z. In Appendix C we provide figures 
similar to Fig. [T] constructed for different colour pairs. 

Photometric measurements in the it band have rela- 
tively poor quality compared to g and r, therefore the resid- 
uals of the (NUV — r,u — r,M r ) relation are about four 
times larger than those of (NUV — r,g — r,M r ). An addi- 
tional factor increasing the scatter is an important difference 
between the extinction vector direction and the blue slope 
of the relation in the (NUV — r,u — r) plane compared to 
(NUV — r,g — r) increasing the scatter at NUV — r < 4 mag. 

The (NUV —i,g — i, Mi) relation has fitting residuals 
about 50 per cent higher than the (NUV — r,g — r, M r ) one, 
although one would expect them to be similar given very 
high quality of i and r band photometry and similar depen- 
dence of these colours on the stellar population evolution. 
We explain this by higher uncertainties of the fc-correction 
computation in the i band connected to a broad range of 
the Ha + [Nil] emission line strength in our galaxies. In our 
sample, except the very low-redshift objects (Z < 0.03), 
the observed-frame i band may be contaminated by the 
Ha + [Nil] emission in a galaxy which may vary a lot from 
object to object. However, during the fc-correction compu- 
tation we rely on some average line streng t h pro vided by 
the template spectra of iBlanton fc Roweisl (J20Q3). There- 
fore, objects with very weak or very strong emission lines 
will have their observed r — i colours redder or bluer than 
what is predicted by the templates and what was included in 
the 2D-polynomial approximations of ^-corrections. Hence, 
the i band fc-correction values may be biased. Even though, 
given a ~50 per cent larger range of g — i colours compared 
to g — r, the relation in this colour space may be used in the 
same way as the one presented in Section 3. 

The (NUV — z, g — z) colour combination provides an- 
other very good alternative to (NUV — r,g — r, M r ), however 
with ~80 per cent larger scatter because of lower photomet- 
ric quality in the z band compared to r in the SDSS. The 
remarkable features of this relation are: (a) a notably higher 
luminosity tilt in the red sequence region and (b) low resid- 
uals in the blue part of the relation caused by even better 
coincidence of the colour change direction due to the stel- 
lar population evolution and internal extinction than in the 
(NUV — r, g — r) plane. 

The (NUV — z, r — z) colour pair starts to suffer from 
very similar behaviour of r and z burdened by relatively 
high photometric errors in the z band. Therefore, although 
the r — z colour range is nearly the same as that in g — r, 
the dispersion of the residuals is notably higher that compli- 
cates the usage of the (NUV — z, r — z, M z ) colour-colour- 
magnitude relation. 

4-2.2 Combinations including Far-UV and NIR colours 

GALEX FUV measurements have on average much worse 
quality than the NUV ones because of the lower detector 
sensitivity and also lower fluxes for intermediate-age and old 
galaxies. However, we still can see similar colour-colour- 
magnitude relations if we use FUV magnitudes instead 
of NUV although with higher residuals especially in the 
low-luminosity part of the relation where the mean surface 



brightness of galaxies decreases. The computed fc-corrections 
also have higher uncertainties in FUV as well as internal ex- 
tinction effects introducing additional scatter. 

The red sequence region in the colour-colour projec- 
tion extends from ~2.5 mag in NUV — r to ~4 mag in 
FUV — r due to even higher sensitivity of FUV colours to 
small fractions of young stars. However, all combinations 
involving FUV magnitudes are sensitive to the UV upturn 
in old early-type galaxies (|Codel 11969^ iBertola et al. 1982) 
likely caused by the stellar evolution (|Yi et al.lll997h T which 
results in the ambiguity of the relation in the red sequence 
region. That is, after some "turning point" (7-8 Gyr), the 
FUV — r colour becomes bluer when the stars are getting 
older. 

We used NIR UKIDSS photometry to test the exis- 
tence of photometric relations in the combinations involv- 
ing optical-NIR colours. None of the combinations except 
(NUV — Y, g — Y, My) provides a relation having simi- 
lar tightness to what we detected in the optical colours: the 
fitting residuals are of a n order of 0.2 m ag or larger. 

It is known (see e.g. lMaraston|[2005l ') that JHK colours 
are sensitive to AGB stars presenting in intermediate-age 
stellar populations, and that at certain ages (1-2 Gyr) the 
optical-NIR colours (g — H or r — H) are dominated by 
them being redder than the colours of old stellar popula- 
tions by a few tenths of a magnitude. Then, given a much 
larger range of e.g. NUV — H than that of r — H , this ex- 
cess will be significant, and it will strongly depend on the 
SFH of a given galaxy, so that galaxies with intermediate 
NUV — H colours having different SFH families may have 
significantly different g — H or r — H colours smearing out 
the intermediate-to-red part of the relation except its very 
red end. In addition, optical-NIR colours are more sensitive 
to the metallicity than the optical ones. Hence, the natural 
relatively low metallicity spread of galaxies at a given lumi- 
nosity will introduce high scatter of their g — H or r — H 
colours. Because of the AGB phase, for truncated SFHs, the 
colour evolution in the (NUV — H , g — H) plane will also be 
more complex than in the optical colour and it will strongly 
depend on the truncation time. 



4.3 Redshifts from three photometric points 

This section of the paper aims at an independent mathemat- 
ical proof of the existence of the tight UV-optical colour- 
colour-magnitude relation of galaxies. The detailed discus- 
sion of the technique and its practial applications to the 
existing photometric survey data will be provided in a sep- 
arate paper. 

The mathematical consequence of the relation and 
smooth dependencies of fc-corrections on observed colours 
is the possibility of the existence of a univocal functional 
dependence of a redshift on observed colours and magni- 
tudes of galaxies. Such a dependence, if found, would con- 
firm the existence of the universal colour-colour-magnitudc 
relation. Importantly, it arises from a non-zero curvature of 
the colour-colour-magnitude surface and significantly differ- 
ent colour-magnitude distributions for the two colours used. 
In a degenerated case, e.g. (u — r, g — r, M r ) where the two 
colour-magnitude distributions are very close to the linear 
dependence and, therefore, galaxies reside on a surface very 
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similar to a plane, the photometric redshift determination 
becomes impossible. 

For the following computations we define the two sub- 
samples from the main galaxy sample extended to the red- 
shift z = 0.52 (270 016 objects) by using their fc-corrected 
g — r colours, red galaxies (g — r >0.73 mag. 77 070 objects) 
and blue galaxies (g — r <0.7 mag, 167 157 objects) exclud- 
ing a small fraction of objects in the "green valley". These 
samples were again separated into a low- (z ^ 0.25) and 
high-redshift (z > 0.25) parts containing 214 770 (32 317 
red + 160 310 blue) and 56 275 (45 365 red + 7 227 blue) 
galaxies correspondingly. Here, most of blue galaxies in the 
high-reds hift sample come from the dee p SDSS Stripe 82 
imaging (jAdelman- McCarthy et all 120061 ') and the fraction 
between blue and red galaxies clearly demonstrates the tar- 
get selection algorithm of SDSS biased towards luminous 
red galaxies at intermediate and high redshifts. 

We approximated the spectroscopic redshifts of galaxies 
in our low-redshift subsample as a three-dimensional poly- 
nomial function of observed r, NUV — r, and g — r, and at- 
tempted to recover the photometric redshifts Zphot (Fig(SJ). 
The dispersion a(A(z)) of the residuals A(z) = (z pno ^ — 
-2spec)/(l+zspec) of 0.025 together with catastrophic failure 
rate (defined as fraction of objects with A(z) > 3a(A(z))) 
of <t] = 0.8 per cent is comparable to the best available 
photometric redshift techniques exploiting multi-band FUV- 
to-NIR photometr y, sophisticated mathem atical and sta- 
tistical algorithms (|Wav fc Sri vastava 2006) and additional 
morphological information (jWrav fc Gunn 120081 ) . For red 
and blue galaxies, the residuals and the catastrophic fail- 
ure rates were: a(A(z)) re d = 0.021, rj rc d = 2.2 per cent and 
a(A(z)) b iue = 0.024, ry b i uo = 0.7 per cent. 

Consequently, at higher redshifts when the restframe 
NUV photometric band shifts to the optical domain, it 
should be possible to determine photometric redshifts pre- 
cisely using it — r — z, g — r — z, or g — i — Y broadband pho- 
tometry. We tested this hypothesis with our high-redshift 
subsample by fitting their redshifts as a function of observed 
SDSS (z, u — z, r — z) and obtained the residuals having a 
dispersion a(A(z)) — 0.036 and the rate of catastrophic fail- 
ures rj = 1.1 per cent. For red and blue galaxies, the residuals 
and the catastrophic failure rates were: cr(A(z)) ro d = 0.034, 
?7rod = 1.0 per cent and cr(A(z))biue = 0.047, ?7biuo = 
1.2 per cent. These relatively large errors are mostly due 
to the very poor quality of it-band Petrosian magnitudes 
having typical uncertainties of an order of 0.3 mag for high- 
redshift galaxies. We notice here, that if one uses model 
magnitudes instead of Petrosian ones, the relation becomes 
much tighter for the red subsample (er(A(z)) re d = 0.027), 
however it nearly disappears for blue galaxies whose it-band 
model magnitudes do not correspond to their real photo- 
metric properties because of light distributions being very 
far from regular exponential or de Vaucouleurs profiles. 

The demonstrated possibility of the precise photometric 
redshift computation for both, red and blue galaxies with 
a small fraction of outliers from three photometric points 
involving a near-UV and optical colours proves the existence 
of the tight unversal colour-colour-magnitude relation for 
normal galaxies of all types, not only red sequence objects. 

We compare these metrics to the existing photomet- 
ric redshifts techniques. Using the "le phare" photometric 
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Figure 5. Recovery of galaxy distances from three photometric 
points. Density histograms of photometric redshift {Zp^ ^) de- 
termination by fitting their spectroscopic redshifts (Zspec) as a 
polynomial function of three parameters: an observed magnitude 
and two observed colours corrected for the Galactic extinction. 
The blue plot displays a sample of 214 000 low-redshift galaxies 
(0.03 < Z < 0.25) with NUV, g, and r photometric measure- 
ments from GALEX and SDSS. The green plot is sample of 56 000 
intermediate-redshift galaxies (0.2 < Z < 0.52) whose redshifts 



J phot 



were determined using u, r, and z SDSS DR7 photometric 



data. Red and blue solid lines denote median of Z^oi and the 
standard deviations of {Z^ot ~ ^spec)/(l + ^spec) residuals for 
blue and red galaxy subsamples respectively. 



redshift code combined with a template optimisation proce- 
dure and the application of a Bayesian approach, based on 
the sample of galaxies with 9 individua l photometri c mea - 
surements of a quality similar to ours, lllbert et al.1 (|2006l ) 
find the dispersion a(A(z)) to be 0.025 and r\ = 1.9 per cent 
(though we stress that the authors defined the catastrophic 
failure limit at the fixe d level of 0.15 wh i ch cor responded to 
6cr in their statistics). iMobasher et al.l (|2007l ) analyse the 
performance of different photometric redshift codes on a 
dataset that comprises 16 photometric points for every SED 
in question. The best result achieved in their study with 
own method reaches the dispersion of residuals as low as 
a(A(z)) = 0.033 and r\ = 2.2 per cent. The SDSS database 
provides sev eral photometric redshift e s timat es obtained as 
described in lAdelman- McCarthy et al.l (|2007h . We have ex- 
tracted three of them: (1) those obtained from the compar- 
ison of the observed colors of galaxies to a semi-empirical 
reference set (hereafter photoz) from the photoz table, and 
(2) neural network estimators derived from galaxy magni- 
tudes ("Dl") and (3) colours ("CC2") from the photoz2 
table. In general, they perform quite well for red galaxies 
in our high-redshift subsample (cr(A(z)) ro d, photoz = 0.020 
and 77 r od, photoz = 1.7 per cent, o-(A(z)) rcd ,Di = 0.020 
and ijrcd.Di = 1-6 per cent, cr(A(z)) red ,cc2 = 0.023 and 
?)rod,cc2 = 1.7 per cent). However, similarly to our approach, 
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the quality is worse for blue galaxies: o-(A(z)) b iuc, P hotoz = 
0.042 and ?7biuc,photoz = 3.4 per cent, o-(A(z))bi U e,Di = 0.039 
and r?biuo,Di = 2.0 per cent, cr(A(z))bi uc ,cc2 = 0.044 and 
»7CC2 = 1.7 per cent. One has to keep in mind that (1) these 
techinques use a lot of additional information (e.g. morphol- 
ogy, size, etc.) and (2) all the galaxies in the spectroscopic 
SDSS sample in fact constitute a training sample of these 
methods, so one has to check those numbers against massive 
third-party redshift surveys. 

We also note that since the redshift determination by 
three photometric points is a mere mathematical conse- 
quence from the colour-colour-magnitude relation, photo- 
metric redshift outliers are themselves objects that fall aside 
from the relation, namely PSGs, dusty starbursts, and AGN. 

To summarise, compared to existing photometric red- 
shift techniques, presented method requires a factor of 3 to 
5 more modest investment in observing time (due to the 
fact that individual galaxy measurements do not need to be 
made in as many photometric bands), being able to provide 
redshifts for large samples of galaxies at the same or bet- 
ter level of accuracy. Moreover, proposed polynomial evalu- 
ation is significantly simpler from the methodological point 
of view than \ 2 minimization with Bayesian priors used in 
mainstream photometric redshift codes. 

There are two main disadvantages of our approach: (1) 
it is not precise for non-typical galaxies, i.e. outliers from 
the colour-colour-magnitude relation; (2) it works only for 
those regions of the parameter space, that are well sampled 
with spectral redshift measurements. Latter means e. g. that 
it is possible to go beyond the SDSS spectral sample magni- 
tude limit retaining the declared precision of our method if 
one uses an external source of spectral redshifts to calibrate 
the functional relation to work with the SDSS photometry 
at fainter magnitudes. But without such a calibration, pho- 
tometric redshifts estimates for faint galaxies will be wrong. 

Conceptually, the presented multi-dimensional polyno- 
mial fit resembles the training of artificial neural networks 
som etimes used for the ph otometric redshift determination 
(e.g. iD'Abrusco et al.ll2007T ) , though the underlying machin- 
ery is different. In both cases, there is a non-linear transfor- 
mation (a 3D-polynomial function in our case or consequent 
multi-level sigmoid transformations in case of neural net- 
works) of some input measurements into the output redshift 
estimate. And the coefficients of a transformation are tuned 
( "trained" ) in a way to work as good as possible for the ref- 
erence ("training") dataset. Hence, both methods unlike \ 2 
template fitting family are insensitive to systematic errors 
of the input data. That is, the "templates" are constructed 
from the dataset itself and are tolerant to its problems by 
construction. 



5 SUMMARY 

We presented the universal very tight colour-colour- 
magnitude relation in optical and near-UV filters followed 
by the vast majority of non-active galaxies of all morpho- 
logical types covering at least 8 magnitudes in luminosity 
from the sample including 225 000 low-redshift (Z < 0.27) 
galaxies observed by SDSS and GALEX surveys. A spe- 
cial case is the connection of the optical g — r colour to 
the NUV — r colour and M r luminosity which we approxi- 



mated by a low-order polynomial surface with the residuals 
of ~0. 03-0. 07 mag in the entire covered luminosity (—23.5 < 
M r < -15.5 mag) and colour (0 < NUV - r < 7.5 mag) 
ranges. 

We have demonstrated that there is a strong correlation 
between the NUV — r colour and the galaxy morphology 
while for the optical g — r colour this correlation is much 
weaker. We identified several classes of the outliers consti- 
tuting about 2 per cent of the total galaxy population and 
explained their nature, the most important being rare com- 
pact elliptical and poststarburst E+A galaxies. 

We have performed stellar population modelling and 
shown that the relation can be explained by galaxies having 
constant and exponentially declining SFHs, while truncated 
SFHs with very rapid colour evolution after the SF cessation 
explain the properties of E+A outliers. The most important 
conclusion is the predicted weak dependence of the colour- 
colour-magnitude relation shape on a redshift up-to Z ~ 0.9 
that will allow to use it for search of distant galaxy clusters 
in 3-colour broad band images. 

The existence of such a tight photometric relation sug- 
gests the possibility of the high-precision photometric red- 
shift estimates using only three photometric points that we 
have also illustrated. Our empirical photometric redshift 
technique using a minimalistic input dataset has a quality 
comparable to or better than most published Z^ Q ^ tech- 
niques, while it is much simpler to implement. 

The detailed astrophysical interpretation of the de- 
scribed photometric relation still has to be done. However, it 
can be used already as a powerful galaxy classification and 
selection instrument based only on their integrated photom- 
etry as well as a tool to search for representatives of rare 
galaxy types in photometric galaxy samples. 
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APPENDIX A: VALIDATION OF THE RESULT 

Two factors may in principle lead to the spurious creation 
of the colour-colour-magnitude relation presented in this 
work: (a) sample selection biased toward specific colours cor- 
responding to the described surface, (b) serious faults in the 
fc-correction computation artificially bringing most galaxies 
to that relation. 

As far as we apply no selection to the SDSS DR7 
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spectroscopic sample of galaxies based on their morphol- 
ogy, colours, sizes etc., and GALEX is a full-sky survey 
providing photometric information for all detected objects, 
the only source of selection effects may be the spectro- 
scopic target selection of SDSS. T he corresponding algo- 
rithms are exhaust ively described |Eisenstein et al] l200ll ; 
IStrauss et ail l2002h . and we find no evidences for them 
to introduce any biases on galaxy selection that may lead 
to the creation of "empty" regions in the colour-colour- 
magnitude parameter space. This is also well illustrated 
by the broad and very well filled colour distribution of 
SPSS galaxies (IStrateva et all l200ll ; iBlanton et ail l2003al ; 
iBaldrv et al.ll2004 ) a nd studies of galaxy l uminosity func- 
tions based on SDSS (|Blanton et al.ll2003bl '). 

In order to test the existence of potential problems in 
the fc-correction computations, we have conducted two spe- 
cific tests. For this, we selected two sub-samples of galaxies 
in narrow redshift ranges, Az%: 0.03 < z\ < 0.05 (24 319 
galaxies) and Az 2 : 0.08 < z 2 < 0.10 (37303 galaxies). We 
chose relatively low-redshift samples because the SDSS tar- 
geting algorithm has a limiting magnitude r — 17.77 mag in 
a 3 arcsec aperture so that galaxies at higher redshifts do 
not sample well the luminosity axis of the parameter space. 

The first test was to compare the colour-colour- 
magnitude distributions of galaxies in the two sub-samples. 
In case of fc-correction computation problems, one would 
expect systematic differences between the two sub-samples, 
which we have not detected (see Fig. lAlfa for an "edge-on" 
view of the relation) . 

The second test was to entirely disable the fc- 
correction computation for these two sub-samples. Because 
fc-corrections do not change significantly within the nar- 
row redshift sub-samples, but do change between the sub- 
samples, one would expect to get two tight colour-colour- 
magnitude sequences qualitatively resembling the relation 
for fc-corrected magnitudes, but with some quantitative dif- 
ferences such as the blue cloud slopes and red sequence 
positions. We obtained the result exactly as predicted by 
this intuitive assumption. The edge-on views of the non fc- 
corrected relations are presented in Fig. lAlb . 

However, the strongest argument supporting the ex- 
istence of a universal colour-colour-magnitude relation in 
UV-optical colours is the possibility of computing photo- 
metric redshifts using only three observed colours as demon- 
strated in the Discussion section. 



APPENDIX B: FITTING SURFACES INTO 
STRONGLY NON-UNIFORM 
THREE-DIMENSIONAL SCATTERED 
DATASETS 

A surface fitting procedure is an essential mathematical 
component required to obtain results presented in the paper. 
The observational photometric datasets for galaxies from 
wide-field survey have strongly non-uniform distribution in 
the colour-colour-magnitude space due to the superposi- 
tion of the complex distribution of galaxies connected to 
their physics and various selection effects and observational 
biases. 

Using visual inspection, we revealed a distribution of 
points in the colour-colour-magnitude space close to a 
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Figure Al. Tests with the two sub-samples of galaxies in narrow 
redshift ranges. The top panel demonstrates the number density 
plots for galaxies in two narrow redshift ranges in the projection 
shown in green and blue, corresponding to the edge-on view of 
the colour-colour-magnitudc relation. The bottom panel is the 
same as the top one, but without applying fc-corrections to the 
measurements. 



smooth surface, but applying standard x 2 -b ase d linear sur- 
face fitting techniques did not yield the results of a reason- 
able quality because: (a) the density of points in the NUV- 
colour-magnitude plane varies by several orders of magni- 
tude while the individual measurements have comparable 
quality; (b) distribution of points around the surface is some- 
times significantly asymmetric and non-Gaussian. The for- 
mer property of the distribution leads to the fact that the 
scarcely populated regions can deviate significantly from the 
surface without notable change of the goodness-of-fit as the 
best-fitting surface tries to minimize the deviation in the 
densest regions of the parameter space. The latter property 
leads to the biased fitting results as the x 2 technique as- 
sumes the Gaussian distribution. The K-sigma clipping tech- 
nique will not solve the problem here because we are dealing 
with a large number of points deviating from the symmetric 
distribution and not with individual outliers. 
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To tackle these issues in a simple way, we decided to 
use a two-step technique for the surface fitting. First, we 
defined a fine grid (cell size of about 0.25x0.25 mag) in the 
NUV-colour-magnitude plane and in every bin computed 
median values of the optical colour being fitted. This al- 
lowed us to pick up the maxima of the (e.g. g — r) colour 
distributions in every bin and not the mean values, that was 
critical in order to account for the asymmetrical distribution 
of points around the surface. Then we filtered out the val- 
ues where the 2D-histogram counts in a given bin were below 
some threshold (usually, 5 or 7 galaxies). At the second step, 
we fitted a low-order polynomial surface into these median 
values using a standard routine fitting linearly the polyno- 
mial coefficients and assigning equal weights to all points 
remained after the filtering at the first stage. This way we 
took into account the strongly non-uniform distribution of 
galaxies in the NUV-colour-magnitude plane. 

This two-step approach resulted in almost flat distri- 
bution of residuals displayed in Fig. [T] and in the next Ap- 
pendix. 



APPENDIX C: 

COLOUR-COLOUR-MAGNITUDE RELATIONS 
IN DIFFERENT COLOUR COMBINATIONS 

In Fig. IC1HC4I we show the colour-colour projections of 
the colour-colour-magnitude relation in different near-UV- 
optical colours described in Section 4.2. All the fitting results 
including the coefficients of the best-fitting surfaces, fitting 
residuals in colour-magnitude bins and other essential in- 
formation for the usage of the relations are provided in the 
electronic forrrQ. 
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Figure CI. The same as Fig. [T] but for the (NUV — r, u — r, 
M r ) space. 



4 http:/ /specphot. sai.msu.ru/galaxies/ 
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Figure C2. The same as Fig.[T]but for the (NUV -i, g-i, Mi 
space. 



NUV-z, mag 

) Figure C3. The same as Fig. [T] but for the (NUV - g, g 

M z ) space. 
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Figure C4. The same as Fig. [T] but for the (NUV — z, r — z, 
Mz) space. 
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